Images of the Dark Soliton in a Depleted Condensate 
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The dark soliton created in a Bose-Einstein condensate becomes grey in course of time evolution 
because its notch fills up with depleted atoms. This is the result of quantum mechanical calculations 
which describes output of many experimental repetitions of creation of the stationary soliton, and 
its time evolution terminated by a destructive density measurement. However, such a description 
is not suitable to predict the outcome of a single realization of the experiment were two extreme 
scenarios and many combinations thereof are possible: one will see (1) a displaced dark soliton 
without any atoms in the notch, but with a randomly displaced position, or (2) a grey soUton with 
a fixed position, but a random number of atoms filling its notch. In either case the average over 
many realizations will reproduce the mentioned quantum mechanical result. In this paper we use 
A'^-particle wavefunctions, which follow from the number-conserving Bogoliubov theory, to settle 
this issue. 
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I. INTRODUCTION 

Let us start with a trivial but important remark on 
measurement in quantum mechanics: in general a single 
measurement on an iV-particle system is neither equiva- 
lent to many measurements on the iV-particle system, nor 
to N measurements on N one-particle systems. The only 
exception is the case of a product iV-particle state where 
all particles are described by the same single-particle 
wavefunction. Not surprisingly the particles must be 
bosons forming a Bose-Einstein condensate, and a mea- 
surement on this A'^-particle state can be replaced with 
N measurements on N one-particle systems in the same 
single-particle wavefunction. 

The importance of this " nonequivalence" principle can 
be seen in the explanation of an interference experiment 
Ref. in an excellent paper by Javanainen and Yoo , 
and then elaborated in Refs.|3l|. 
consider a Pock state 



In Ref.l2i the authors 



N N 
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with N/2 atoms in a condensate of wavefunction (f>Q{x) 
and N/2 in a condensate with an orthogonal wavefunc- 
tion (f>i{x). A single particle density matrix of the state 
([T|) predicts that on average a density measurement will 
give a distribution of atoms 



(2) 



without any interference fringes between the two con- 
densates. This single particle distribution is an average 
over many density measurements on the same Pock state 
. Simulation of a destructive density measurement 
based on the full TV-particle wavefunction and the 



actual experiment , give distributions with an interfer- 
ence term 



N 



( |0o(x)|2-)-|0i(x)|2 + 

+MxWi{x)e+'^ + c.c.) 



(3) 



For a given experimental realization the phase ip is cho- 
sen randomly from the interval [0, 27r). The outcome of 
a density measurement is the same as if before the mea- 
surement all N atoms were not in the state but in an 
absurd single condensate with a wavefunction 



(4) 



with a relative phase ip that is not known before we ac- 
tually measure the density distribution. It turns out Q 
that it is enough to measure only a small fraction of the 
large total N to prepare the remaining atoms in one of 
the condensate wavefunctions or, in other words, to 
establish a definite phase ip. Destructive measurement of 
atomic positions, which effectively annihilates measured 
atoms from the trap, drives the remaining atoms into 
a condensate state with a randomly picked condensate 
wavefunction. These interesting conclusions cannot be 
obtained without full knowledge of the A^-particlc state 
(^: the single particle density distribution is a (mis- 
leading) average over many experiments with different 
outcomes ip. To predict possible outcomes of individual 
measurements on an A'^-particle state it is essential to 
know its TV-particle wavefunction. 

A perfect condensate is a state where all iV atoms are 
in the same single particle state described by a single 
particle wavefunction (j)o(x), 



(5) 



When we include interactions between atoms but neglect 
depletion of atoms from the condensate wavefunction, 
then 00 is a solution of the celebrated Gross-Pitaevskii 



2 



equation Including small quantum depletion in the 
framework of the number-conserving Bogoliubov theory 
leads to pair- correlated eigenstates . For exam- 
ple, the state without any quasiparticles (the Bogoliubov 
vacuum) has a pair-correlated form 



(dt)f |0) 



|0) 



(6) 



Operators ak annihilate in a basis of functions (pk orthog- 
onal to the condensate wavefunction (j)Q. This general 
pair-correlated ansatz © has been known for a while, 
see e.g. the review by Leggett 0, but a general solution 
for the coefficients and corresponding wavefunctions 
4>k has been found only very recently in Ref.j^] (the so- 
lution in the special case of a uniform condensate was 
given in Ref.Q). The two-particle creation operator d'' 
has to commute with all quasiparticle annihilation oper- 
ators (see Appendix A) 



= 0, for aU TO. 



(7) 



The TV-particle state (O is a foundation on which one 
can build theory of BEG entirely in the language of iV- 
particle wavefunctions. 

In our two recent papers 0| we studied depletion 
from a condensate with a dark soliton. The dark soliton 
state is a collectively excited condensate with a notch in 
the condensate wavefunction. In the Thomas-Fermi limit 
of strong interactions the notch behaves like a soliton 
with a wavefunction close to the notch proportional to 



(poix) ~ tanh(a; — X) , 



(8) 



where X is a position of the soliton measured in units 
of the healing length. The anomalous Bogoliubov mode 
with a wavefunction 



01 (x) ~ —dx4>oix) = cosh '^{x — X) 



(9) 



dominates quantum depletion from the condensate 
within the soliton. In Ref.j^ we found an average num- 
ber of atoms dN depleted into the non-condensate mode 
and a single particle density distribution which close 
to the soliton a,t X — looks like 

Pi{x) = iN-dN)\Mx)\^+dN\Mx)\'' ■ (10) 

For dN = the distribution has a hole a,t X = 0, but 
with increasing quantum depletion dN the hole is filling 
with atoms. This single particle result means that on 
average over many experiments the notch of a depleted 
dark soliton appears filled with atoms (it appears grey, 
not dark), but it gives no clue what to expect as an out- 
come of a single destructive density measurement. We 
address this issue in the following sections. 

The paper is organized as follows. In Sections |lIIIIIIIVl 
we apply the number conserving Bogoliubov theory 



(shortly summarized in Appendix A) to the dark soli- 
ton, and simulate density measurement on different N- 
particle quantum states that all may pass under the same 
label: "depleted dark soliton" . We find that images of 
the "depleted dark soliton" depend qualitatively on the 
actual A^-particle state. To simulate the measurements 
we adapted the numerical algorithm of Ref.*?] which is 
described in Appendix B. We conclude in Section Ivl 



II. 



STATIONARY DARK SOLITON 



The dark soliton state 10] of a condensate in a ID har- 
monic trap can be defined as an antisymmetric solution 
of the stationary Gross-Pitaevskii equation 

- ^9^00 + ^a;^(?!)o +g|0op0o = /i0o (H) 

In this paper we use the dimensionless oscillator units. In 
the Thomas-Fermi (TF) limit oi g ^ 1 the dark soliton 
state is 



(j)Q{x) w F{x) tanh(a;/Zo). 



with 



Fix) 



l2n 
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(12) 



(13) 



where the chemical potential fi ~ {3g/2)'^^^ /2 and the 
healing length ~ V2{2/'igf'^. 

A perfect condensate with all TV atoms in the dark soli- 
ton mode 00 is not an energy eigenstate because inter- 
actions between atoms continuously deplete atoms from 
the condensate wavefunction 00 • The Bogoliubov theory 
is a systematic way to describe small quantum fluctua- 
tions around the Gross-Pitaevskii solution, and to find 
stationary states of a depleted condensate. 

Solution of the Bogoliubov-de Gennes (BdG) equa- 
tions (see Appendix A) reveals many different modes of 
the quasiparticle spectrum. However, only the negative 
energy (anomalous) mode ® has a wavefunction local- 
ized in the soliton notch 0, |3, 0| ■ iFrom all modes the 
anomalous mode will contribute the most to the density 
of incoherent atoms in the soliton notch. In the follow- 
ing, similarly as in Ref. "gll, we truncate the Bogoliubov 
spectrum to the anomalous mode alone and look for the 
Bogoliubov vacuum state in the form 



|0h:iV)-(44-f Aa^aD* |0) 



(14) 



with a\ creating in the mode 0i. The anomalous mode 
solution of the BdG equations H35() with i 



IS 



ui = - (/+ + /-) , v^ = - iU - /-) , (15) 



with functions 
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VScosh^ (f 



Fix) , (16) 
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see Ref.|2l- For g ^ 1 equations H15|l are dominated by 
/+, (/_|_|/+) (/-I/-), and the two functions become 
the same, ui w ui w As a result in the TF hniit 

the anomalous mode does not mix with other modes and 
the dominant eigenfunction is 



(17) 



as has been verified numerically in Ref. 6]. To find the 
asymptotic behavior of A for g 3> 1 we have to be more 
careful and still keep /_ for a moment. The eigenvalue 
can be obtained from the general prescription presented 
in Appendix A and it is 



A 



(/+!/+) -(/+!/-) 

(/+!/+) + (/+!/-) ^ g'/' 



(18) 



Knowing A we can find an average number of depleted 
atoms dN in the Bogoliubov vacuum state H14|l . In the 
particle Fock representation the unnormalized vacuum 
state H14II is a sum 



N) ^ ^A' 



fc=0 



V(iV- 2fc)! (2fc)! 
(iV/2- fc)! kl 



N~2k,2k)(19) 



For depletion greater than unity but small with respect 
to the total number of atoms N, i.e. 1 ^ dN <C N, we 
may use the Stirling formula, and then replace the sum 
by an integral. As a result we find expressions for dN 
and its dispersion D{dN) 



dN'. 



1 ^ g'/' 
h-iX^ A 2V6 31/3 
D{dN) = dN. 



(20) 
(21) 



The single particle density distribution in the vacuum 
state (|14I19|I is 



Pi{x) ^{N- dN)\Mx)\^ + dN\Mx)? 



(22) 



The first term on the RHS is the dark soliton density 
profile with a hole at a; = 0, while the second one stands 
for the depleted atoms localized in the soliton notch. A 
good measure of depletion is a ratio of atomic density in 
the notch at a: = to the density near the notch 



dN\<f>iiO)\' 
NF{0)^ 



N 



1/3 



(23) 



suggesting that for sufhciently large number of atoms, 
the notch will not be visible. For the parameters of the 
Hannover experiment T^], where g ~ 7500 and N ~ 
1.5 X 10^, dN « 60 atoms and the ratio is 15%. 

However, Eq. (|22|l predicts an average over many ex- 
perimental realizations. For a single experiment one has 
to choose randomly positions of N atoms according to 
the A'^-particle probability distribution. In general this 
is rather a formidable task, but, as shown by Javanainen 



and Yoo [2j (see Appendix B), for a state spanned on only 
two modes, and with a sequential algorithm that chooses 
a position of n+ 1-st atom after positions of n atoms have 
already been "measured" , this problem is polynomial in 
N. 

In general we expect that a single measurement of a 
depleted dark soliton will show the same atomic density 
distribution as if all N atoms were initially prepared in 
a condensate with one of wavefunctions 



(j)o{x) - (j)i{x)qe'"^ 



(24) 



The real parameters {q, cp) are chosen at random and fluc- 
tuate from experiment to experiment. In the special case 
of (/3 = 0, TT the condensate wavefunction H24|l is simply a 

displaced dark soliton 



4>o{x) — (f>i{x)q w F{x) tanh 



y V 4 



(25) 



In general, when ^ 0, tt, the condensate wavefunction 
(|24|l gives a density distribution 

N\^o-Me"'\'' = 

N [00 - (fiiq cos (pf + Ncflq^ sit? . (26) 
For small q this distribution can be interpreted as a soli- 
ton shifted to q^ ^ cos (/s, but with a nonzero atomic 

density N(j)\q'^ siv? Lp in the notch. For 93 7^ 0, tt the soli- 
ton appears as a displaced grey soliton. 




FIG. 1: Density measurement for state |06 : A''), Ea. l|19|l . The 
histogram represents number of particles in bins around the 
solitonic notch measured in a single experiment. The solid line 
stands for p{x) = A\(j)o{x)\'^ + B\(l}i{x)\'^ + C(l}o{x)(l}i{x), where 
A,B,C are parameters of ISUII averaged over all atoms, see 
the formula 15UI I. The dashed line gives the expected average 
outcome of many experiments i.e. the single particle density. 
The values of the parameters correspond to the Hannover 
experiment [ij, where g = 7500 and iV = 1.5 x 10^. 

Adopting the scheme of Javanainen and Yoo we gen- 
erated results of individual density measurements of the 
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Bogoliubov vacuum state H14I19() . A generic outcome is a 
displaced dark soliton shown in Fig.Q] In Fig.|21we show 
how definite values of cos{(p) and q are established in the 
course of measurement of subsequent atomic positions. 
It is enough to measure only a small fraction of atoms to 
suppress fluctuations of ip and q, see Fig. El 
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FIG. 2: Values of cos Lp (upper plot) and q during the density 
measurement performed on the state |0b : A*'). 

A density measurement of the Bogoliubov vacuum 
state (I14ll9f) always gives a displaced dark soliton with 
(/3 = in Ea. H24|l . At first sight this may seem obvious 
because the wavefunctions 0o and 0i are real, and all 
amplitudes in the state H19|l are also real. However, this 
simple argument does not take into account the possi- 
bility that the measurement of a small fraction of atoms 
may localize the state of remaining atoms not in a single 
condensate (|24|) . but in a superposition of several conden- 
sates ;3| with different parameters {q,(p). Indeed, as we 
can see in Ea. H26|) . density distributions for {q,+ip) and 
{q, —If) are the same so the density measurement cannot 
distinguish ~\-ip from —tp.K real superposition of two con- 
densates with opposite phases has real amplitudes, just 
like the state (|19|l . We checked that a measurement of a 
Fock state \N — n^n) , instead of the Bogoliubov vacuum 
state (|14ll9f) . indeed results in a real superposition of two 
condensates with nonzero +Lp and —95, which appears as 
a displaced grey soliton. What makes a difference is the 
fact the Bogoliubov vacuum (|19|l is a superposition over 
many Fock states. It is the phase coherence between dif- 
ferent Fock states, all with real amplitudes, that enforces 
the observed ip — Q. 



III. SOLITON CONDENSATE WITHOUT 
INITIAL DEPLETION 

In current experiments 0, the dark solitons are 
generated with the help of the phase imprinting . As 
we argued in Ref. 0, when the depletion in the conden- 
sate before the imprinting is negligible, then the soliton 
state right after imprinting can be idealized by the state 



without any atoms depleted from the solitonic conden- 
sate 00 to the anomalous mode 



\^{t = 0)) = \N,Q) . 



(27) 



This state is not a stationary state. To obtain time evo- 
lution of this state we consequently truncate the Bogoli- 
ubov Hamiltonian to the anomalous mode where 



(28) 



with 



hi = 



V^(l-A2) 



(29) 



With this Hamiltonian the initial perfect condensate 
\^[t = 0)) evolves into a depleted state \^{t)). In Fig. 01 
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FIG. 3: Single particle densities of the state |*(t)), Eq.lgJJ, 
at different instants of time. The values of the parameters 
correspond to the Hannover experiment (T^. where g « 7500 
and iV = 1.5 X 10^ 



we show how a single particle density evolves in time. 
After about 15 ms (for the parameters of the Hannover 
experiment 12J) the hole in the single particle density 
distribution fills up with depleted atoms. 

We simulated single realizations of destructive density 
measurements of the state \'^{t)). A generic density dis- 
tribution after evolution for 10 ms is shown in Fig. 0] 
We see that the soliton is dark and it is displaced with 
respect to the trap center. Inspection of the state after 
10 ms shows that {N-k, A;|*(10 ms)) ~ e''"^ with a tiny 
if « 0.01. In principle the time evolution might generate 
any phase (f € [0,27r), but the parameters of the exper- 
iment ^3 such that there is not enough time to get 
any more significant ip from the idealized initial state H27|) 
before the linearized Bogoliubov theory breaks down. 
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FIG. 4: The same as FigClbut for the state = Q.15)>. 
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FIG. 5: The same as FigQbut for the state |rand). 



IV. GENERIC NON-STATIONARY SOLITON 
STATES 

In a realistic experiment the dark soliton is created in 
neither the stationary nor the zero depletion state dis- 
cussed in the previous two Sections. The two states may 
appear "natural" to a theorist's eye accustomed to such 
highly idealized examples. However, these states are not 
realistic from experimental point of view. The soliton 
is phase imprinted on a condensate with a 10% thermal 
cloud i.e. in a highly excited non-stationary state. We do 
not know what the actual state, or ensemble of states, is, 
but we can still test how robust is the picture developed 
in the previous Sections with respect to imperfections of 
the quantum state. 

To get an idea what are images of a generic non- 
stationary state we model the state by a state with ran- 
dom phases. The Hamiltonian of the system is diagonal 
in the Fock basis of quasiparticles so we expect a generic 
non-stationary state in to have apparently random phases 
in particle Fock representation. To be more specific, we 
take as an example the Bogoliubov vacuum H19|l with 
randomized phases a^. 



I rand) 



N/2 
k=0 



^ {N - 2k)\{2k)\ 
{N/2-k)\k\ 



\N-2k,2k) . (30) 



Single particle density in this state is exactly the same 
as for the Bogoliubov vacuum, see Eq. (|22|) , but possible 
outcomes of density measurements can be very different 
in the two states. Indeed, a generic outcome in the state 
|rand) shown in Fig. ]E\ reveals a displaced grey soliton 
with cos(^ = 0.45. The phase disorder present in the 
state |rand) translates into random tp's. As a second ex- 
ample, we take the state |rand) as an initial state instead 
of the state (|27|l . evolve it for 0.7 ms, and then measure 
atomic positions. In Fig|Slwe show a single particle den- 
sity after 0.7 ms and a corresponding generic outcome of 
density measurement. Not surprisingly, we again see a 



displaced grey soliton, this time with cos ip — — 0.- 
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FIG. 6: The same as Fig^but for the state |rand) after evo- 
lution over 0.7 ms. 



These two examples demonstrate that phase disorder 
in particle Fock representation can change qualitatively 
possible results of density measurements. 



V. CONCLUSION 

In Refs.[^|3 we have considered influence of quantum 
fluctuations on a density profile of a Bose-Einstein con- 
densate excited to a solitonic state. We have shown that 
the soliton notch may be filled with atoms that are de- 
pleted from the condensate due to interactions between 
particles. The analysis has been performed by calculat- 
ing a single particle density matrix. However, such a 
matrix provides informations about the density profile 
that is a result of average over many experimental real- 
izations. To obtain predictions for a single experimental 
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realization one has to choose randomly positions of atoms 
according to the multiparticle probability density. 

In this paper we concentrate on the analysis of out- 
come of the destructive density measurement in a sin- 
gle experiment. Applying the number conserving Bogoli- 
ubov theory of a BEC we calculate TV-particle quantum 
states of a depleted condensate with a dark soliton. The 
iV-particle wavefunctions are used to see what happens 
after the depleted condensate is subject to an ideal de- 
structive (density) measurement of atomic positions. We 
have found that: 

• For a BEC excited to a stationary state, i.e. the 
Bogoliubov vacuum state without any quasiparti- 
cles, a destructive measurement reveals density dis- 
tribution which appears like a displaced dark soli- 
ton with an empty notch. Average displacement is 
proportional to the average number of incoherent 
atoms in the notch. 

• For the non-stationary state that evolves from an 
initial state without any incoherent atoms in the 
soliton notch (i.e. a perfect condensate with a dark 
soliton) a destructive measurement reveals, simi- 
larly like in the previous case, a displaced dark 
soliton. The longer the evolution lasts the larger 
the depletion of the condensate becomes (and thus 
also average displacement of the soliton). 

• The previous two cases are related to highly ideal- 
ized states of a BEC. To analyze influence of imper- 
fections (in the preparation of the soliton conden- 
sate) on the result of density measurement we in- 
troduced random phases in amplitudes of different 
Fock states corresponding to the Bogoliubov vac- 
uum state. In such a non-stationary state the mea- 
surement reveals a displaced grey soliton. Both 
the average displacement and the average density 
of atoms in the notch of this grey soliton are pro- 
portional to the number of atoms depleted from the 
condensate. 

Very different A^-particle quantum states may corre- 
spond to the same single particle density matrix. How- 
ever, different states result in qualitatively different re- 
sults of destructive density measurements. In view of the 
strong dependence of measurement results on a quantum 
state, it is an open question what A^-particle state is rel- 
evant for the present day dark soliton experiments |l2j| . 
The phase imprinting method ,141] prepares neither an 
energy eigenstate nor the initial state without any de- 
pleted atoms. Certainly, further experiments would be 
helpful to settle this problem. 

Quantum fluctuations of a condensate are an intrigu- 
ing result of interatomic interactions that introduce two- 
particle and higher correlations into the state of the sys- 
tem. BEC as a correlated system might be a very in- 
teresting toy model for condensed matter physics. De- 
pletion to the notch of the dark soliton gives a unique 



opportunity to detect quantum fluctuations simply by 
density measurement. A dark soliton in a strictly one- 
dimensional trap is well within the reach of present tech- 
nology. 

APPENDIX 

A. The N-particle State of BEC 

Here we briefly summarize main conclusions of the 
number-conserving Bogoliubov theory of BEC A 
field operator is split into a condensate part aQ(j)Q{x) 
and a non-condensate part 5tp{x), 

■ip(x) = 4>o{x)an + (5-0(f) . (31) 

To zero order in the "small" fluctuation operator 5i}j{x) 
the A^-particlc state is a perfect condensate 

{alY |0) (32) 

with all N atoms in the condensate wavefunction (j)o{x). 
The (po solves a stationary Gross-Pitaevskii equation 
(CPE) 

= -^V20o + V{x)(l)o + g|<^^|0o • (33) 

To zero order the interaction affects only the shape of 
the condensate wavefunction 0o (x) through the nonlinear 
term ^I^qI^o in the GPE. 

Quantum depiction from this condensate shows up in 
the second order of the perturbation theory where the A^- 
particle state becomes a pair-correlated Bogoliubov vac- 
uum state 

N_ 

\Ot:N)^{d^)^ |0) = ialal+J2^kalal\ |0). 

(34) 

The sum runs over an orthonormal basis of non- 
condensate modes 4>k{x) orthogonal to (j)f){x). The eigen- 
values Afc and the eigenmodes (j)k {x) have been calculated 
in full generality only very recently 0. Their construc- 
tion runs in the following steps. Its full justification can 
be found in Refs.0,0|. 

• Solve the stationary GPE to get 4>o{x) and ^. 
Note that (/)o does not need to be the ground state, 
it can be an excited stationary state with a dark 
soliton or vortex. 

• Solve the linear Bogoliubov-de Gennes equations 
for the Bogoliubov normal modes of small fluctua- 
tions around 00 j 

-\v^Um + V{x)U^ + 2g|</.opC/™ + g(t>lV^ = 



7 



(35) 



to get the eigenvalues ujm and "raw" Bogoliubov 
eigenmodes Um{x), Vm{x). 

• Project the raw [/,„, Vm on the subspace orthogonal 

to 00, 



Vr, 



= U„i - 4'a{4'o\U,n) , 

= Vrn - 4'n{4>o\Vm) , 



(36) 



to get the Bogoliubov eigenmodes {u„i,Vm). 
• Normalize the modes {u„i,v,-a) so that a norm 

(Wm|Um) - (Wm|Wrn) = +1 • (37) 

Ignore modes with negative norm. 

The normalized Bogoliubov modes define bosonic 
quasiparticle annihilation operators 



t 



(38) 



which transfer atoms between the condensate and the 
non-condensate modes, but conserve the total number of 
atoms N . For small depletion the Hamiltonian can be 
approximated by a sum of harmonic oscillators 



H = 



E 



(39) 



The A^-particle Bogoliubov vacuum H34|l is the eigenstate 
of H without any quasiparticles, 



b,^\%:N)=Q 



(40) 



The A^-particle state (|34|) is annihilated by all 6's when 
the defined in Ea. H34(l commutes with all 5's, 



= 



(41) 



In general construction of such a d) runs along the fol- 
lowing steps. 

• Choose an orthonormal basis of states <j>k{x) or- 
thogonal to 00 J and then calculate matrices U and 



V with elements 

U„ik = {Um\4>k) , 
Vmk = (0fe|Wm) ■ 

• Invert U and then calculate a matrix 
Z = IJ-^ V . 



(42) 
(43) 



(44) 



• Diagonalize Z to get its eigenvalues Afc and or- 
thonormal eigenmodes (pu ■ In the basis of (pk the Z 
becomes a Z = diag{Ai, A2, . . .}. 

Having calculated Afc and 0^ we can construct the oper- 
ator and the Bogoliubov vacuum state (|34() . Excited 
eigenstates can be obtained by repeated action of the 
quasiparticle creation operators on the vacuum state. 



B. Density Measurement 

In an ideal density measurement positions of all N 
atoms are measured. The probability distribution for the 
N positions in a A'^-particle quantum state \^) is given 
by 



Pn{xi, ■ ■ ■,xn) ^ 

m'^^Xi) . . .i;\xN)i^{xN) 



.ij{xi)\^) (45) 



We want to simulate density measurements by generating 
from this distribution typical outcomes {xi, . . . ,X]sf)- In 
general this is a formidable task. 

The problem becomes tractable when, for some phys- 
ical reasons, it is possible to truncate the single par- 
ticle Hilbert space to just two modes, say, (f)o{x) and 
<j>i{x) (generalization to 2,3... modes is straightfor- 
ward). With only two modes we can adopt the algorithm 
of Javanainen and Yoo 0. In their algorithm positions 
{xi, . . . ,xn) are not generated all at once but one after 
another. To begin with, xi is chosen randomly with a 
reduced single particle probability distribution 

Piixi) ^ m\xi)il>{xi)\-^) . (46) 
Once the actual xi is chosen we calculate a state vector 



l^-i) = V(a;i)|^') = [aoM^i) + ai(f>iixi)] l*> 



(47) 



and store it in the memory. After (k — 1) steps we know 
(xi, . . . , a;fc_i) and a state vector 



l*fc 



k-l 



tp{xk-l)\'^k-2) 



(48) 



calculated for the actual Xk-i- To generate the next co- 
ordinate Xk we use a conditional probability distribution 



p{xk I Xk-l 



Pk{Xk,- . . ,Xi) 



Pk-l{xk~l, 

^k-MKxk)Mxk)\^k-i) 
(*fc_ll*fc_l) 



,Xl) 



(49) 



Note that this is a function of Xk only, all coordinates 
Xk-i, ■ ■ ■ ,xi have already been fixed. Once we get the 
actual outcome Xk we calculate \^k) = '^Axk)\^ k-i) and 
store it in the memory. 

The key simplification discovered in Ref.0] is that 
all single particle probability distributions that we en- 
counter, like Eas. (|46l49|) . have the same functional form 

p[x) = A \Mx)?+B \Mx)?+ 
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C cj)*{x)M^) + C* M^Wii^) ■ (50) 

To find the real A, B and the complex C we need to 
sample p{x) at only four values of x. When (j)o,(j)i are 
real, like for the dark soliton, then C is also real, and 
just 3 sampling points are enough. The total algorithm 
to generate N positions (xi, . . . , x^v) scales like N'^, a 
slight improvement over the in Ref.|3]- 

To summarize, we provide an input A^-particle state 
I (in a Fock representation build on the two single par- 
ticle modes (/)o, </>!) and the algorithm generates a random 
string of coordinates (xi,...,xjv) according to the A^- 
particle probability distribution H45fl . Each application 



of the algorithm returns a random outcome {xi, . . . ,Xk) 
of a single realization of density measurement on the state 

l*>- 
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